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Słowa kluczowe 
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Streszczenie 


W pracy przedstawiono koncepcję probabilistycznego modelowania wystą- 
pienia różnorodnych zagrożeń stanu zdatności systemu technicznego. Głównym 
celem tego opracowania są konstrukcje funkcji służących do analizowania kon- 
kurujących zagrożeń. Praca ta została zainspirowana książką [4], w której roz- 
winięte zostały metody probabilistyczne mające szczególne zastosowania w bio- 
logii i medycynie. W punkcie 1 przedstawione są funkcyjne charakterystyki 
niezawodnościowe zaadaptowane na potrzeby badania niezawodności, bezpie- 
czeństwa i zagrożeń systemów technicznych. Pokazane są związki między nimi. 
Najważniejsze wyniki są zestawione w punkcie 2, gdzie przedstawiono model 
konkurujących zagrożeń i jego funkcyjne charakterystyki. W punkcie 3, podane 
są dwa proste przykłady, w których zastosowano przedstawioną metodę mode- 
lowania probabilistycznego. 


Wprowadzenie 


Eksploatację systemu technicznego uważa się za podstawowy etap spraw- 
dzania jego przydatności i spełniania oczekiwań społecznych. Współczesne 
systemy techniczne coraz częściej są monitorowane w czasie eksploatacji, 
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w celu wykrycia zagrożeń stanu zdatności. Interdyscyplinarne powiązanie pro- 
blemów eksploatacji systemu wyraźnie wskazuje na dominującą wśród nich rolę 
diagnostyki zarówno technicznej, jak i użytkowej. Przedstawione w [1] probabi- 
listyczne metody prognozowania przyczyny utraty zdatności systemu dychoto- 
micznego zostały w tej pracy przeniesione na grunt badania różnorodnych za- 
grożeń zdatności lub trwałości systemu. 

Przyjmujemy, że T jest nieujemną zmienną losową o wartościach rzeczywi- 
stych określoną na jednorodnej populacji systemów 42, tj. T: 2 — R. Jeżeli dla 
danego systemu obserwujemy kilka zdarzeń, to stosujemy wektor losowy. 
W teorii niezawodności zmienna losowa T jest czasem zdatności lub trwałości 
obiektu technicznego, tj. losowym czasem opisującym zdolność systemu do 
realizacji nakładanych zadań, w określonych warunkach eksploatacji. Podsta- 
wową miarą utraty zdatności lub pojawienia się zagrożenia w przedziale czasu 
[0, £] jest powszechnie stosowana w probabilistyce dystrybuanta, tj. funkcja F: 
R-[0, 1], zwana w teorii niezawodności funkcją zawodności, przy czym dzie- 
dzina funkcji zawodności ograniczona jest do R,. Symbol R oznacza zbiór liczb 
rzeczywistych, a R, oznacza zbiór liczb rzeczywistych nieujemnych. 

Aby być zgodnym z normami ISO, przyjmujemy, że dystrybuanta jest funk- 
cją co najmniej prawostromnie ciągłą, tj. F(t) = P(T < t), gdzie P jest miarą pro- 
babilistyczną indukowaną przez zmienną losową 7. 


1. Charakterystyki funkcyjne czasu utraty zdatności systemu 


Niech T oznacza czas do wystąpienia negatywnego zdarzenia w kontekście 
pewnego systemu. Zdarzeniem tym może być trwała bądź chwilowa utrata zdat- 
ności całego systemu lub jego części, pojawienie się lub rozwój jednego lub 
kilku z możliwych zagrożeń itp. Dualnie można wyróżnić pozytywne zdarzenie, 
takie jak ustąpienie zagrożenia, w wyniku podjętych pewnych działań, rozpo- 
znanie zagrożenia, przywrócenie zdatności systemu itp. 

Oprócz dystrybuanty F pięć innych funkcji znanych z teorii niezawodności 
charakteryzuje rozkład zm. 1. 7, są nimi: 

a) funkcja przetrwania S (the survival function), zwana też niezawodnością, 

b) funkcja ryzyka h (the risk function), interpretowana tutaj jako intensywność 
utraty zdatności, 

c) funkcja prawdopodobieństwa lub gęstości prawdopodobieństwa f, 

d) skumulowana funkcja ryzyka H, 

e) oczekiwany pozostały czas zdatności mrl (the mean residual life). 

Funkcja niezawodności S: R,—>[0, 1] jest podstawową miarą bezpieczeń- 
stwa, tj. prawdopodobieństwa niewystąpienia negatywnego zdarzenia lub zjawi- 
ska w przedziale czasu [0, t]. Funkcja ta określona jest wzorem S(t) = P(T > t). 
Jeżeli funkcja S jest absolutnie ciągła, to można ją przedstawić w postaci cał- 
kowej 
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S(t)=] fwdu, 120 (1) 


Funkcja f, spełniająca równanie (1), to gęstość prawdopodobieństwa. Jeżeli 
istnieje gęstość, to można określić funkcję intensywności utraty zdatności. 
W teorii niezawodności funkcja ta jest także znana jako warunkowa intensyw- 
ność uszkodzenia, w demografii znana jest jako współczynnik umieralności, 
w procesach stochastycznych jako funkcja intensywności, a w ekonomii jako 
odwrócony iloraz Milla. Intensywność utraty zdatności jest zdefiniowana 
wzorem: 


POST <t+At|T2t) 


2 
iF (2) 


h(t) = lim 
At>0 


Jeżeli czas T jest typu ciągłego, to h(t) = f(1)/S(t) = —d ln(S(t))/ dt. 


W kontekście konkurujących zagrożeń stanu zdatności, do ustalania przy- 
czyny utraty zdatności systemu służą przyczynowo wyspecyfikowane intensyw- 
ności zagrożeń. Ze wzoru Taylora mamy przybliżenie 


S(t) — S(t + At) 


f@= re 


(3) 


oraz 


S(t)-S(t+ At) 


h = 
4 S(t)At 


(4) 


Ze wzoru (3) wynika, ze gęstość prawdopodobieństwa można interpreto- 
wać jako spadek niezawodności lub bezpieczeństwa w małym przedziale czaso- 
wym o długości Af, a z (4) wynika, że intensywność utraty zdatności to wzgled- 
ny spadek niezawodności w takim przedziale. W populacji obiektów na tyle 
licznej, aby można założyć „funkcjonowanie prawa wielkich liczb”, ze wzoru 
(3) wynika, jaka część obiektów utraci zdatność lub będzie zagrożona w czasie 
At, gdy jako podstawę przyjmiemy liczbę wszystkich zdatnych obiektów na 
początku (t = 0), natomiast za pomocą (4) można określić tę część, w stosunku 
do obiektów zdatnych w chwili t. 

Z funkcją intensywności utraty zdatności związana jest skumulowana funk- 
cja intensywności H, zwana w teorii niezawodności funkcją wiodącą, którą 
można interpretować jako informującą o wyczerpywaniu się „zapasu zdolności 
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do przebywania w stanie zdatności”. Każdą z wymienionych pięciu funkcji cha- 
rakteryzujących czas do wystąpienia zagrożenia obiektu można wyrazić przez 
dowolną inną spośród nich. Zależności te są podane w [2]. 

Jeżeli T' jest zmienną losową typu dyskretnego, to związki pomiędzy funk- 
cjami są inaczej opisane. Przypuśćmy, ze T przyjmuje wartości ¢;, j = 1, 2,... 
z prawdopodobieństwami określonymi przez funkcję prawdopodobieństwa 
ft) = PUT = t), j = 1, 2,..., gdzie tı < to <.... Funkcja przeżycia dla zmiennej 
losowej typu dyskretnego jest określona wzorem: 


St) =P >) = DG) (5) 


tj >t 
Funkcja intensywności utraty zdatności jest określona wzorem 


JC) 
Sj) 


h(t,) = PC =1, {0 >t,)= „J=l, 2,..., gdzie S(t) = 1 (6) 
Ponieważ f(t) = S(G-1) — S(t), więc h(t) = 1 — S)/S(G-1), j = 1, 2,... 


Zauważmy, że funkcję niezawodności można zapisać jako iloczyn warun- 
kowych prawdopodobieństw przeżycia 


SO = JISE) Sta) (7) 


tj St 


Stąd wynika związek pomiędzy funkcją niezawodności a funkcją intensyw- 
ności 


SQ) = [[Q-h(t,)) (8) 


t,śt 


Skumulowaną funkcję intensywności utraty zdatności dla czasu typu dys- 
kretnego można zdefiniować na dwa sposoby. Jeśli przyjmiemy definicję 


H(t) = Dh) (9) 


4 
tj St 


to znany z przypadku ciągłego czasu zdatności związek S(t) = exp{—H(a)} prze- 
stanie być prawdziwy. Aby zachować prawdziwość tego związku, należy przy- 
jac definicję z [3] 
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H=- nd- h(t,)) (10) 


t,śt 


Jeżeli wielkości h(t,) są małe, to (9) można zastosować do aproksymacji 
(10). W zastosowaniach zalecane jest stosowanie definicji (9) z powodu możli- 
wości użycia bezpośredniej estymacji z próby losowej dla danych z ograniczo- 
nym czasem monitorowania. Ponadto otrzymane estymatory mają pożądane 
własności statystyczne. 

Piątą charakterystyką funkcyjną, jaką można zastosować w analizie zagro- 
żeń stanu zdatności, jest pozostały średni czas do wystąpienia zagrożenia 
w chwili ¢, w teorii niezawodności zwany oczekiwanym pozostałym czasem 
przetrwania (mean residual life — mrl). Dla obiektów w wieku ¢ funkcja mri(t) 
jest miarą charakteryzującą ich oczekiwany pozostały czas do wystąpienia za- 
grożenia w chwili £. Funkcja mrl zdefiniowana jest wzorem 


mri(t) = E(T- t|T> t) (11) 


Geometrycznie mrl(t) jest równa polu pod krzywą funkcji przetrwania S(t) 
na prawo od £, podzielonemu przez S(t). Zauważmy, że oczekiwany czas na za- 
grożenie m = E(T) jest równy polu pod krzywą niezawodności, tj. m = mrl(0). 

Dla zmiennej losowej 7 typu ciągłego 


| @-nf(@dx f Sd 


I(t) = 12 
mrl(t) Só 5 (12) 
a dla zm. 1. T typu dyskretnego 
tia OSH) + Dd, 5.4, E54 ESC; 
mrl(t) = aS OEE Byala) eee (13) 


S(t) 


2. Model rywalizujacych zagrozen 


W praktyce często trudno jest wskazać elementarne zdarzenia, które są rze- 
czywistą przyczyną utraty zdatności przez system. Zadanie staje się łatwiejsze 
jeśli ustalimy możliwe źródła utraty zdatności. Na przykład system o strukturze 
szeregowej złożony z k elementów może utracić zdatność na co najmniej k róż- 
nych sposobów. Inżynierowie rozpatrują różne możliwe sposoby utraty zdatno- 
ści jako formy uszkodzenia i często one są rozważane jako pierwotne przyczyny 
uszkodzenia. To ujęcie inżynierskie jest w tej pracy rozwijane. 

Zakładamy, że przyczyną wystąpienia zagrożenia lub utraty zdatności sys- 
temu jest wystąpienie jednego spośród wielu wykluczających się zdarzeń. Mo- 
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del, w którym uwzględnia się wystąpienie pierwszej przyczyny utraty zdatności 
systemu, spośród skończonej ich liczby, nazywamy modelem rywalizujących 
zagrożeń. Model ten znany jest w naukach biomedycznych pod angielską nazwą 
competing risks model, patrz [4]. 

Niech E = (e, é2,..., ex} będzie zbiorem zagrożeń stanu zdatności systemu. 
Przykładem zagrożenia stanu zdatności systemu może być awaria jego elementu, 
błąd operatora, zła konserwacja, nieplanowane obciążenia lub jakiś czynnik 
zewnętrzny jak uderzenie pioruna lub ptaka w obiekt. Do dyskusji nad proble- 
mem rywalizujących zagrożeń formułujemy model z ukrytymi czasami ich wy- 


stąpienia. 
W odróżnieniu od stosowanych metod w inżynierii niezawodności, tutaj 
przyjmijmy, że T; i= 1,..., k jest potencjalnym nieobserwowalnym czasem do 


wystąpienia i-tego konkurującego zagrożenia. Nieobserwowalność zagrożeń 
wynika z ich skrytości działania. Obserwowany jest natomiast stan zdatności 
systemu. W szczególności rejestrowany jest moment wystąpienia zagrożenia 
utraty zdatności systemu i przyczyna tego zagrożenia. Monitorujemy więc czas 
T = min{T; i= 1,..., k} oraz wskaźnik 6 o wartościach w zbiorze {1,..., k}, 
który informuje nas, które z k zagrożeń jest przyczyną utraty zdatności obiektu, 
tj. ô= i jeśli T = T,. Metoda opisana tutaj opiera się na intensywnościach przejść 
ze stanu zdatności do jednego spośród k identyfikowalnych i rozróżnialnych 
stanów zagrożeń zdatności systemu. Podstawową charakterystyką funkcyjną 
rywalizujących zagrożeń jest przyczynowo wyspecyfikowana intensywność 
utraty zdatności systemu ze względu na i-te zagrożenie, zdefiniowana wzorem: 


P(t<T,<t+At,d=i T, 2t, j=1,...,k) 


dla i€ {1,..., k} (14 
AG ie + (14) 


h,(t) = lim 
At>0 


Funkcja Aft) jest intensywnością zagrożenia stanu zdatności systemu 
z i-tego powodu przy założeniu, że system jest narażony na wszystkie zagroże- 
nia. 

Całkowita intensywność zagrożenia stanu zdatności systemu, przy założe- 
niu (2) jest addytywna, tzn. jest sumą wszystkich k intensywności zagrożeń, tj. 


k 
hy (t) = 2 h,(t) (15) 
i=l 


Wyspecyfikowaną intensywność zagrożenia można otrzymać z łącznej 
funkcji bezpieczeństwa dla k rywalizujących zagrożeń, definiowanej wzorem: 


S(t,..., ty) = P(T1> ty,..., TĘ> ty) (16) 


4-2009 PROBLEMY EKSPLOATACJI 13 


Mianowicie 


SQ Gs SE fo ee (17) 
SGD 


h,(t)= 


A przy założeniu, że zagrożenia zdatności systemu występują niezależnie 


— 08, (£,)/0t, 
O) 


h,(t) = dla t,=t (18) 


Oznacza to, że dla niezależnych rywalizujących zagrożeń brzegowe inten- 
sywności zagrożeń i przyczynowo wyspecyfikowane intensywności zagrożeń są 
identyczne. W przypadku zależnych zagrożeń zasada ta nie obowiązuje. Poka- 
zane to jest na przykładach w punkcie 3. 

W modelu konkurujących zagrożeń może wystąpić potrzeba przyjęcia 
pewnych założeń dotyczących współzależności pojawienia się potencjalnych 
zagrożeń. Jeżeli przyjmiemy, że obserwujemy jedynie czas do wystąpienia za- 
grożenia systemu oraz jego przyczynę, a nie czasy wystąpienia potencjalnych 
zagrożeń, to niestety założenia o współzależności zagrożeń nie będą weryfiko- 
walne. Problem ten nazywa się dylematem identyfikowalności współzależności. 
Dylemat ten wynika stąd, że mając dane typu (T, 6), nigdy nie rozpoznamy pary 
zależnych zagrożeń z pary niezależnych obserwacji. 

W problemach konkurujących zagrożeń, obok intensywności zagrożenia, 
może interesować nas prawdopodobieństwo uwzględniające naszą wiedzę 
o możliwości wystąpienia tylko pewnej grupy zagrożeń, a nie wszystkich. 
W tym celu wprowadzamy prawdopodobieństwa dla trzech różnych informacji 
o zagrożeniach. Nazwiemy je prawdopodobieństwem: brutto, netto i częściowo 
brutto. Prawdopodobieństwem brutto nazywamy prawdopodobieństwo zagroże- 
nia stanu zdatności systemu ze szczególnej przyczyny, w rzeczywistym świecie, 
gdzie wszystkie zagrożenia oddziałują na system. Prawdopodobieństwem netto 
nazywamy prawdopodobieństwo stanu zdatności systemu, w hipotetycznym 
świecie, gdzie specyficzne zagrożenie jest jedynym zagrożeniem oddziałującym 
na system. Jest to prawdopodobieństwo brzegowe dla specyficznego zagrożenia. 
Prawdopodobieństwem częściowo brutto nazywamy prawdopodobieństwo utra- 
ty zdatności systemu, w hipotetycznym Świecie, gdzie część zagrożeń została 
wyeliminowana. 

Prawdopodobieństwo brutto jest wyrażane poprzez funkcję skumulowanego 
specyficznego zagrożenia, zdefiniowaną wzorem F;(t) = P(T < t, = i). Funkcję 
tę można wyznaczyć bezpośrednio z funkcji łącznego skumulowanego zagroże- 
nia, łącznej gęstości zagrożeń lub z intensywności zagrożenia specyficznej przy- 
czyny. Otrzymujemy więc zależności 
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FQ) = PO t,6 =i) = [hywexp-H,(u) du (19) 
0 


Tutaj H(t) jest skumulowaną intensywnością zagrożenia stanu zdatności 
systemu, tj. 


kt 
H,(t)= > [h,(u)du (20) 


J=10 


Zauważmy, że funkcja F(t), ie {1,..., k} zależy od intensywności, w któ- 
rej występowały wszystkie k konkurujące zagrożenia, a nie tylko i-te. Ponieważ 
intensywność zagrożenia h{t) może być estymowana bezpośrednio z obserwo- 
wanych danych, więc zawodność F((t) jest również bezpośrednio estymowana, 
bez czynienia żadnych założeń dotyczących łącznego rozkładu wystąpienia za- 
grożeń. Funkcja F;(t) nie jest dystrybuantą, ponieważ F;(ce) = P(Ó= i). Funkcja 
ta jest niemalejąca oraz F;(0) = 0 i Fico) < 1, więc nazywamy ją semi- 
dystrybuantą. 

Funkcja netto bezpieczeństwa S;(ż) jest brzegową funkcją prawdopodobień- 
stwa, wyznaczoną z łącznej funkcji bezpieczeństwa, przyjmując, że t;= 0 dla j #i. 
Jeżeli rywalizujące zagrożenia są niezależne, to funkcja bezpieczeństwa netto 
jest związana z prawdopodobieństwami brutto poprzez formułę 


+ dF (u) 
S.(t)= -| —— d 21 
KO) eaf a u (21) 


Jeżeli zagrożenia są zależne, to prawdopodobieństwa bezpieczeństwa netto 
można oszacować poprzez prawdopodobieństwa brutto, tj. 


Sp) <S;(1) $1- F(t) (22) 


W przypadku prawdopodobieństw częściowo brutto przyjmujemy, że J jest 
podzbiorem pozostałych potencjalnych zagrożeń utraty zdatności przez obiekt, 
natomiast J© jest podzbiorem zagrożeń wyeliminowanych. Niech T” = min{X,, 
ie J}. Definiujemy semi-dystrybuantę częściowo brutto 


F’! =P(T! <t, =i), ieJ (23) 
Tutaj i-ta semi-dystrybuanta przedstawia prawdopodobieństwo utraty zdat- 


ności obiektu z powodu i-tego zagrożenia, jeśli obiekt jest narażony jedynie na 
zagrożenia ze zbioru J. 
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Na zakończenie zdefiniowana jest intensywność częściowo brutto zagroże- 
nia zdatności 


— 9S (t-t) Ot, 
S (tistg) 


4 (1) = oraz t; = t, gdy j EJit,=0, gdy j EJ” (24) 


Stąd semi-dystrybuanta częściowo brutto przedstawia się wzorem 


F? ()=P(T" <t,=i)= jal (x) a- >[2 add |as (25) 
0 


jeJ 0 


Jeżeli zagrożenia stanu zdatności są niezależne, to intensywność zagrożenia 
częściowo brutto upraszcza się do postaci 


dF.(t)/ dt 


y= 
F S;(t) 


(26) 


3. Przykłady 


Przykład 1. Przypuśćmy, że stan zdatności obiektu jest zagrożony ze 
względu na trzy niezależne przyczyny, które mogą wystąpić zgodnie z rozkła- 
dem wykładniczym o intensywnościach /,, A2, 43, odpowiednio. W tym przy- 
padku intensywności wystąpienia zagrożenia netto i brutto są takie same. Z ad- 
dytywności, łączna intensywność utraty zdatności obiektu hr = A; + A, + As. 
Funkcja zagrożenia zdatności obiektu z powodu i-tej przyczyny ma postać 


4, 


TARŁA 


(1-exp-L(4 +4, + A;))) (27) 


Zauważmy, że prawdopodobieństwo brutto utraty zdatności systemu 
z pierwszej przyczyny, w czasie £ nie jest takie samo jak prawdopodobieństwo 
netto (brzegowe) utraty zdatności w tym czasie, dane przez 1 — exp(—A,f). Za- 
uważmy, że F;(oo) = A/(A, + A + As). Jest to prawdopodobieństwo, że i-te zagro- 
żenie będzie pierwszą przyczyną utraty zdatności. Jeżeli rozważymy hipotetycz- 
ną sytuację, w której występują tylko pierwsze dwa zagrożenia utraty zdatności 
systemu, tj. J = (1, 2), to intensywność zagrożenia częściowo brutto ma postać 
AFO = X, i = 1, 2, natomiast semi-dystrybuanta częściowa brutto 


F; (t) = | A,exp(-u(A, + 4,))du = 4 (1-exp tA +2,)) +12 (28) 
0 A, +A, 
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Przykład 2. Przypuśćmy, że przyczyną utraty stanu zdatności systemu są 
dwa rywalizujące zagrożenia o łącznej funkcji bezpieczeństwa postaci: 


S(t,,t>) =1+ O(At, +41,)]'/,620,4,,4,20 (29) 


W tym przypadku dwie potencjalne chwile wystąpienia zagrożenia są sko- 
relowane. Współczynnik 7 Kendala wynosi 0/(0 + 2). Z (18) dla ty = > = t jest 


Clearer aa 
Ponieważ T = min(7), 7>), więc dla t; = fh = t jest 
S(,1) =S,()=[1+6t(4 +4,)]'/,620,4,4,>0 (31) 
i intensywność zagrożenia systemu wynosi 
A +A, (32) 


1+0t(4,+4,))"'* 


Zauważmy, że brzegowa funkcja bezpieczeństwa dla 7, jest określona wzo- 
rem 


S(,00=+0t4)"" (33) 


Stąd brzegowa intensywność zagrożenia ma postać /4/(1 + 0 441), a funkcja 
zagrożenia stanu zdatności i-tą przyczyną jest postaci 


2, 


F. = 
TES 


aroraa’) (34) 


Podsumowanie 


W pracy uwaga zostata skupiona na charakterystykach funkcyjnych stoso- 
wanych w modelowaniu parametrycznym zagrożeń stanu zdatności systemu. 
Pominięte zostały metody coraz szerzej używane w nieparametrycznym mode- 
lowaniu. W analizie zagrożeń metodami parametrycznymi łatwo jest interpreto- 
wać przedstawione funkcje i ich parametry, w szczególności intensywność za- 
grożenia. Jeśli wystąpienia zagrożeń są niezależne, to przedstawiony model 
konkurujących zagrożeń można zastosować do wszystkich znanych z teorii nie- 
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zawodności rozkładów prawdopodobieństwa. W przypadku zależnych zagrożeń 
trzeba użyć rzadziej stosowane metody probabilistyki wielowymiarowej. Kon- 
strukcje rozkładów wielowymiarowych mogą odegrać kluczową rolę w rozwoju 
zastosowań probabilistycznego modelowania konkurujących zagrożeń. 
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Summary 


This paper describes a concept of probabilistic modelling of different 
competing risks of a technical system ability state. The main purpose is the 
construction of functions to analyse competing risks. In Point 1, functional 
characteristics are introduced, adopted to investigation of reliability, safety and 
the threats of technical systems. Moreover, relationships between them are 
shown. The most important results are shown in Point 2, where the model of 
competing risks is introduced. In Point 3, there are given two simple examples 
of systems. For these systems, there were appointed the hazard rates, the crude 
sub-distribution functions for the particular competing risks, the partial crude 
hazard rates, and the partial crude sub-distribution functions. In this publication, 
attention is concentrated on functional characteristics applied in parametrical 
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modelling of the threats of the system ability state. In analysis of competing 
risks with parametrical methods, the introduced functions and their parameters 
are easily interpret, in particular, the hazard rate. If occurrences of threats are 
independent, then the introduced model of competing risks can be applied to any 
probability distribution. In case of dependent competing risks, we should use 
methods of multidimensional probability methods, which are seldom applied. 
The multidimensional distributions can play key role in development of 
applications of competing risks analysis. The whole of study finishes short 
recapitulation and cited literature. 


